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Abstract 


A mathematical  model  of  a simplified  atmosphere  is  described 
which  is  suitable  for  considering  certain  processes  of  heat  exchange. 
The  model  is  set  up  for  an  almost  incompressible  atmosphere,  that 
is,  the  variations  in  density  are  neglected  except  in  the  term 
representing  gravitational  force.  Only  the  dynamics,  not  the  thermo- 
dynamics, of  the  problem  is  considered.  The  resulting  equations  of 
motion  are  simplified  by  perturbation  methods  regarding  the  ratio 
of  atmospheric  thic’sness  to  the  radius  of  the  earth  as  the  perturb- 
ation parameter.  Given  a density  distribution  in  the  atmosphere  as 
an  arbitrary  function  of  latitude  and  height,  the  model  determines 
the  zonal  velocity  component  and  the  meridional  circulation.  A 
simple  numerical  example  is  given  to  illustrate  the  type  of  flow  to 
be  expected  from  the  model. 
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Illustrations 

1 . Introduction 

Although  much  effort  has  been  devoted  to  constructing 
mathematical  and  mechanical  models  of  the  atmosphere  with  the  aim 
of  exnlaining  the  general  features  of  the  atmospheric  circulation, 
there  is  at  present  no  satisfactory  manner  of  linking  the  observed 
behavior  of  the  atmosphere  to  the  fundamental  sources  of  energy 
responsible  for  that  behavior.  The  ultimate  source  of  energy  for 
atmospheric  circulation  must  bo  solar  radiation.  The  motion  in  the 
atmosphere  derives  its  energy  from  the  unequal  distribution  of  heat 
produced  by  the  different  races  at  which  heat  energy  is  absorbed 
from  and  radiated  into  outer  space  at  different  parts  of  the  atmos- 
phere. The  motion  in  turn  furnishes  processes  of  heat  exchange 
which  tend  to  produce  a uniform  distribution  of  heat. 

The  circulation  is  observed  to  be  near  a quasi-equilibrium 
state  at  all  times  so  that  large  surpluses  or  deficits  of  heat 
energy  in  the  atmosphere  do  not  occur.  There  are,  however,  signif- 
icant departures  of  the  circulation  from  the  mean  state  which  do  not 
seem  to  depend  u^on  a particular  distribution  of  heat  energy  in  the 
atmosphere  but  rather  suggest  a dynamic  instability  of  the  basic 
zonal  circulation  in  at  least  some  of  its  modes. 

The  problem  of  the  stability  of  zonal  flow  is  of  great  im- 
portance and  has  received  a good  deal  of  attention  in  the  past  few 
years.  Kuo  (1950),  for  example,  has  considered  the  stability  problem 
in  considerable  detail  and  has  discussed  a possible  role  which 
stability  plays  in  maintaining  and  regulating  the  general  circulation. 
Studies  of  the  stability  problem,  however,  do  not  give  very  much  in- 
sight into  how  the  general  circulation  is  maintained  since  the 
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conversion  cf  hoot  into  kinetic  and  potential  energy  is  not  con- 
sidered . 


Much  of  the  difficult-/  in  constructing  a theory  of  the 
general  circulation  is  encountered  in  choosing  a simple,  meaningful 
model  in  which  the  thermodynamics  of  the  atmosphere  may  be  Included. 
The  mathematical  models  which  have  been  developed  for  the  study  of 
dynamic  stability  of  zonal  circulation  are  not  suitable  for  the 
thermodynamic  problem  because  most  of  them  either  do  not  explicitly 
include  the  effects  of  vertical  structure  in  the  atmosphere  or  are 
applicable  only  when  the  thermodynamic  processes  are  negligible. 

The  classical  Hauley  model  is  more  suitable  but  is  difficult  to 
construct  mathematically. 

In  the  classical  Hadley  model  it  can  be  shown  by  qualitative 
reasoning  that  at  least  three  meridional  circulation  cells  must 
exist  to  correspond  approximately  to  the  observe*!  zonal  winds.  The 
subtropic  and  arctic  cells  arc  driven  thermally  while  the  middle 
cell  is  driven  indirectly  by  the  other  tv/o  cells  through  lateral 
friction.  The  arguments  leading  to  this  picture  cf  the  circulation 
consisting  of  three  cells  and  the  dynamics  of  the  middle  cell  have 
been  given  by  Rossby.  A semi-empirical  model  of  a three  cell  type 
has  been  constructed  by  Dorodnitsyn,  Izvekov,  and  Schwetz  (1939) 
from  the  observed  temperature  distribution  in  the  atmosphere  and 
the  surface  pressure  distribution,  using  a method  developed  by 
Kochin  (1935).  Rossby  (19^1)  gave  an  excellent  discussion  of  the 
Hadley  model  and  pointed  out  some  of  the  difficulties  encountered 
in  attempting  to  explain  certain  features  of  the  circulation  with 
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this  model 
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A model  of  a simplified  atmosphere  is  constructed  here  in 
order  that  sene  processes  of  heat  exchange  ear.  be  .studied  analytical- 
ly. -‘'he  purpose  of  this  model  is  to  give  the  distribution  of  zonal 
and  meridional  velocity  components  corresponding  tc  a given  density 
field  and  representing  a balance  of  pressure-gradient*  Coriolis, 
gravitational,  and  frictional  forces.  Assuming  the  density  field 
to  be  known,  the  velocity  field  is  determined.  Using  the  energy 
equation  and  the  equation  of  state,  the  distribution  of  heat  sources 
and  sinks  necessary  to  maintain  the  partici  lar  density  field  in  the 
steady  state  can  be  calculated.  "he  thermodynamics  of  the  model 
can  be  considered  as  a separate  problem  and,  therefore,  is  not  taken 
up  in  detail  here. 


2 . Simplifying  Assumptions 
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The  equations  of  mass  and  momentum  conservation  serve  as  the 
basis  for  the  model  considered  in  the  following  analysis.  It  Is 
necessary  to  make  simplifying  assumptions  at  the  outset  in  order 
to  obtain  a model  sufficiently  simple  for  analytical  treatment. 

The  atmosphere  is  considered  to  be  an  almost  incompressible 
fluid  covering  a smooth  earth  to  a mean  depth  h which  is  very  small 
in  comparison  with  the  earth’s  mean  radius  a , The  motion  of  the 
atmosphere  is  considered  in  the  steady  state  and  is  assumed  to  be 
driven  by  a knov/n  density  dis  iribution.  It  is  assumed  that  the 
density  variations  are  small  enough  to  bo  sufficiently  accurately 
represented  in  their  effect  on  the  fluid  motion  by  considering  a 
variable  gravitational  force.  Frictional  forces  are  assumed  to 
arise  from  Reynolds  stresses  and  are  taken  into  account  by  introducing 
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a constant,  isotropic  kinematic  eddy  viscosity.  In  adopting  an 
eddy  viscosity,  it  is  implicit]  y as  sum  i that  t/.c  stoi  d"  .-notions 
arc  arrived  at  bv  ar.  averaging  process.  The  non-linear  terr.s  arc 
neglected  to  simplify  the  mathematical  analysis  so  tent  conclusions 
based  on  this  model  would  bo  valid  only  when  these  terms  are  rela- 
tively snail.  The  motion  is  assumed  to  be  independent  of  longitude 
since  otherwise  the  motion  docs  not  anpear  to  have  any  obvious 
physical  significance  unless  surface  irregularities  of  the  earth  are 
taken  into  account. 

The  coordinates  of  a point  in  the  atmosphere  may  be  specified 
by  giving  the  distance  z or  die  point  from  the  geoid  measured  along 
the  local  vertical  direction,  the  latitude  9,  and  the  longitude  X 
at  the  geoid.  It  is  assumed  that  this  coordinate  system  can  be 
sufficiently  well  approximated  by  a spherical  coordinate  system 
(X,  9,  r)  in  which  gravitational  forces  act  in  the  radial  direction. 
Here,  X is  the  longitude  measured  eastward  from  a reference  meridian, 
9 the  latitude  measured  from  the  equatorial  plane 'and  r measured 
from  the  origin  at  the  center  of  the  earth.  The  radial  distance  r 
of  a point  at  height  z above  the  geoid  is  approximated  by  a + z, 
where  a is  the  mean  radius  of  the  earth.  The  velocity  components 
(u,  v,  w)  are  chosen  sc  that  u is  the  eastward  component,  v the 
northward  comronent,  and  w the  upward  component  of  velocity. 

3 • Equations  _of  J tot  ion 

On  the  basis  of  the  assumptions  above,  the  equations  of 
motion  can  be  written  in  the  following  form 


2Q  [ co j y v;  - sin  9 v ] = K [ y c' u - — 1 (1) 

V'  cos^  9 


?Qsin  9 u - - .A  3J2  j-  K[  n^v 
pr  09 


v ^ .-*  3w  j 

r^iv"‘  9 r~  69 
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p cr 
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2w  2 3v 
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r r 3 9 
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2 tan  cp  v 
r2 


_1 

r cos  9 


3 cos  9 v 

— 59 


+ J*  JilLs 

r^  3 r 


= 0 


(3) 

(4) 


where  & is  the  earth's  angular  velocity,  p is  the  density  and  is 
considered  constant  while  j(9,r  ),  the  unit -mass  gravitational  force, 
is  modified  to  include  the  effects  of  small  variations  of  the  density. 
The  pressure  is  given  by  p,  and  the  ’’inonatic  eddy  viscosity  by  K, 

The  Laplscian  operator  in  tills  coordinate  system  is 


r^cos  9 39 


..  A.  [ cos  9 -A] 
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1 _3.  r 
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A stream  fv  action  4 is  introduced  for  the  meridional  cir- 
culation so  that 


v = 


v = - 


r cos 
1 


1 34- 

os  9 "3r 


r cos  9 Oy 


(5) 

(6) 


The  equations  can  be  expressed  more  compactly  by  introducing 
0 

the  quantities  M,  x,  D , and  0,  where 


M = r cos  9 u 
v = sin  9 
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„2 


= 3-l.  + a 

O T 

3r*“  9x' 


i = xi+  (X^rl  _a 
3r  r 3x 


end  M is  the  unit-mass  relative  angular  momentum  of  the  zonal  flow, 

p 

and  D , G are  commutative  differential  operators.  The  equations 
become 


- 2 QGOJO  = K D20-i)  (7) 

+ 2 fl  y f?  = - C1  - X2,),  3_p  4-  K a [p2,^]  (8) 

r pr  3x  r 3r 

= D2(*)]  (9) 

r p 3r  r2  3x 


Eliminating  the  pressure  from  (8)  and  (9)  yields 


2 QG(M) 


(1 


- x2) 


is 

3x 


+ K 


D4(^) 


(10) 


The  boundary  conditions  at  the  earth's  surface  are  that 
normal  and  tangential  velocity  components  are  zero.  At  the  upper 
surface  the  cornnonents  of  torque  are  zero.  The  condition  of  con- 
tinuity requires  zero  net  volume  transport  in  the  steady  state  a- 
cross  all  latitude  circles.  In  terms  of  and  i\>  , these  boundary 
conditions  become 


Tj,  = = 0 I 

ar  L 

f at  r = a, 


i':  = 0 J 
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?nd 


ib  = A. 
* a r 


1,  All 

rc-  9r 


_a_ 

3r 


j 


= 0 


« 0 


at  r = 


+ h, 


The  mathematical  problem  is  snecified  once  g is  given  as 
a function  of  x and  r. 


4 . Dimensional  Analysis 

Unfortunately,  it  is  difficult  to  solve  equations  (7)  and 
(10)  regardless  of  the  choice  of  g.  Further  simplification  of  the 
equations  is  achieved  by  dimensional  analysis.  The  equations  are 
converted  into  r.on-dimensicnal  form  by  introducing  the  non-dimen- 
sional quantities  aQ,  z1  , M’ , vj- 1 , g1,  and  p’ , where 

"o  " h/a 

z'  = z/h  or  r = a (1  + a0z'  ) 

M'  = M/aU 
\i) 1 = 4>/p 
g*  = g/S 
P'  = p/ps  h 

and  U,  *P,  g are  characteristic  values  of  the  zonal  velocity  component, 
stream  function,  and  gravity  respectively. 

By  choosing 

<P  = K?  /4Q  2 
U = aQ/V  /2Q 

e = K/2Q  h2 
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equations  (7)  and  (10)  can  be  written  in  the  following  non-dimen- 
sional fori",  the  prir.es  now  doing  omitted 


- 

7:  (4*  ) - 7 

(31) 

G(h) 

= (1 

- x2)  |_S  + e2Dlt’(4>) 

(12) 

where  the  differential  operators  in 

non-dimensional 

form  are 

D2  = 

O 

r. 

3. 

+ a2  -11 

2n  a2 

- x .)  d 

3? 
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+ o.Q  z)  dx^ 

n _ -v- 

ji  — X 

3 

3z 

+ ^ {t- 

- x2)  a 
+ a0z)  3x 

The  order  of  mag 

nitude  of  the  various  terms 

can  be  estimated 

by  choosing  representative  values  of 

the  dimensional 

constants.  For 

example,  the  values 

a 

- 

6371 

km 

h 

- 

10 

Ion 

K 

= 

1C? 

2 -1 

cm  sec 

g 

= 

lo3 

-2 

cm  sec 

P 

= 

10  "3 

g cm**3 

a 

= 7.3X10"5 

sec"-^ 

yield  the  following  characteristic  constants 

U 

- 

110 

m sec”^ 

<P 

- 

7xie9 

3 -1 

mJsec 

ao 

- 

0.0016 

e 

= 

C. 0007 
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5.  Perturbation  E.nuatiqns 


The  investigc.tion  of  equations  (11)  and  (11)  can  he  simpli- 
fied considerably  by  regarding  a0  as  a perturbation  p-rar.eter  and 
expanding  the  solutions  and  the  differential  equation  as  perturba- 
tion series  in  o.Q, 

The  functions  M,  and  g are  expanded  as  power  series  in 


o * ao  M1  + Go  ‘ ‘2  + 


+0  + ao  ’•’l  + a0  +2  + 


* So  + co  Ex  + «0  32  + 


These  series  are  substituted  into  equations  (11)  and  (12)  and  all 

terms  am  collected  as  coefficients  of  terms  in  a power  series  in 

aQ.  Considering  aQ  as  an  arbitrary  parameter,  the  power  series  can 

have  a sum  function  identically  zero  only  if  all  the  coefficients 

of  powers  of  aQ  are  zero.  This  condition  leads  to  an  infinite  set 

of  differential  equations  for  II.  and  4>  . In  the  following  analysis 

J J 

only  the  equations  corresponding  to  the  zero  and  first  order  terms 
in  aQ  are  considered. 


The  equations  of  zero-order  in  aQ  are 

e 2 - x fHo  = - ( 1 - x2 ) £!© 

Sz^  9z  5x 


(13) 


dn0 

"a"? 


with  the  boundary  conditions 


(14) 
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i'o  = tyQ/dz  = o 


^ at  z - C, 


•;0  * c ! 


and 


% - = 0 


3M0/az  = o 


I 


at  z = 1, 


The  equations  of  first-order  in  aQ  are 

2 „ . (i  . x2) 


« -r-  - * 

3z^  dz 

32M 


a m,  5\l',  0 cnu 

+ x --1  = - Cl  - x2)  —2 

^7J~  dz  9x 


aMo  + dsl 
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with  the  boundary  conditions 


4^  = av|),/dz  = ol 

1 1 ^ at  z = 0, 

• ‘1  = o J 


and 


= 


o,  a24>1/dz2  = 2 &4>0/3z"| 
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at  z = 1. 


aMj/  3z  = 2MC 


J 


(15) 

(16) 


The  perturbation  equations  for  the  pressure  can  be  obtained 
in  a similar  manner  from  equations  (8)  and  (9)  after  converting  them 
into  non-dimensional  form. 


^2  _ _ p 

3z  " so 
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xM0  + e2 


d3^o 

3z^ 


(17) 


(18) 
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(19) 

(2C) 


6 . Sol  ujtions.  of  _tho_  Zero -or  de  r_  E pujatipns . 

Solutions  to  the  zero-order  equations  can  be  obtained  in 
closed  form  for  the  special  case  that  3g0/Sx  is  a function  of  x 
only;  so  for  the  purposes  of  studying  simple  flov/s  predicted  by  the 
present  model,  this  special  case  is  considered  in  detail,  "he  more 
general  case  where  0go/9x  is  a function  of  both  z and  x can  be 
solved  in  terms  of  an  integral  over  a Green's  function  or  by  bound- 
ary-layer  methods. 

The  solutions  to  the  zero-orde^  equations  may  be  written 
in  the  form 
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k) 


o = - W— -^2  0(2, 

yC  9x 


M « .CL=j3  *2  » (z,k) 

O x ox  ’ 


where 


and  X j ( j 


0(z,k)  = 1 - T.  Ai  exn  (X  z) 
j=l  3 3 

$ ( z,  ?;)  = f (8(  T,k)  dr 
(o 

1,2, 3, '4),  represent  the  four  roots  of  the  equation 
X^  = - x2/e2  = - 4 k*4'. 


(21) 

(22) 

(23) 

(24) 
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The  roots  are 


\ = - k(l  - i) 

X^  =r  + k(  1 + i) 

X2  = - k(l  + i) 

X4  = + k( 1 - i) 

where 

i = (-1)1/2 

The  function  0(z,k)  can  be  found  by  solving  the  four 
simultaneous  algebraic  equations  arising  from  the  boundary  condi 
tions  for  the  coefficients  Ay  The  function  is  found  to  be 


0(z,k)  = 1 


I 

^ 2 
I sinh  "Sk  - sin  2k 


1 

7 


^ [cos  k cosh  k(l-z)  sin  kz 


- cosh  k sinh  kz  cos  k (1-z)] 

+ [cosh  k cos  k - sinh  k sin  k ] [ cosh  k(l-z)  sin  k(l-z)  ] 

- [cosh  k cos  k + sinh  k sin  k 1 [ sinh  k(l-z)  cos  k(l-z)  ] 

(25) 


There  are  two  cases  for  which  this  expression  is  closely  approximated 
by  a much  simpler  form: 

Case  I.  k > 5 

0(z,k)~l  - (2)^  exp(-kz)  cos  (kz  - nA)  - exp(-k(l-z))  cos  k(l-z) 


0(  z1  k)  - ,(^z)2(l.-zJ.(3“2z) 


Case  II.  k<l 
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Frcfiles  of  <D(z.k)  and  0(z,k)  rre  given  for  several  values  of  k in 
fig.  1 and  fig.  ? respectively. 

The  zero-order  pressure  distribution  can  be  calculated  from 
equations  (17)  and  (lC).  If  9g0/9x  is  a function  of  x only, 
gQ(z,x)  nay  bo  written  as 

g0(z,x)  = gol(z)  + g02 (x) 

and  from  equation  (17),  the  pressure  pQ,  as 

PD(z,x)  = Pol(z)  + go2(*)  C1  + *l(x)  - z ] (26) 

The  equation  for  q(x)  is  obtained  by  substituting  from  (21),  (22) 
and  (26),  the  expressions  for  4»0,  MQ,  and  po  respectively,  into 
equation  (18).  The  resulting  equation  is 


f • 

f 

_ — 

-j 

1 

®S02I» 

88o2 

cosh  2k  + cos  ?k  - 2 cosh  k cos  k 

1 

f 8x 

Ox 

± 

k(sinh  2k  - sin  2k)  _ 

(27) 


The  boundary  condition  for  n(x)  is 

1 


v 


q(x)  dx  = 0, 


The  function  q(x)  represents  the  departure  of  the  upper  surface 
from  the  surface  z = 1 necessary  to  bring  the  pressure-gradient 
forces  into  equilibrium  with  Coriolis  and  frictional  forces. 

The  horizontal  pressure  gradient  can  be  expressed  os  follows 
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'i 


ST  * [ 802(1  + n - z)l 
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= i!o2 

dx 


cosh  2k  4 cos  2k 
k vsi  nh  2"k 


2 cosh  k cos  k 
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It  is  clear  that  the  horizontal  comroncnt  of  the  pressure  gradient 
must  change  sign  in  the  lower  frictional  layer  and  will  be  zero 
along  the  surface 

z = cosh  ?k  + cos  2k  - 2 cosh  k cos  k (29) 

k (sinh  2k  - sin  2k) 

for  any  choice  of  9gQ/3x  which  depends  only  on  x.  Along  this  sur- 
face the  Coriolis  and  frictional  forces  are  of  equal  magnitude. 

7.  Surface  Torque 

It  is  of  interest  to  note  that  equation  (1*0  can  be  integrated 
with  respect  to  z directly.  The  resulting  equation  is 

8M 

Sz * + = 0 <30> 

This  implies  that  the  meridional  transport  at  a given  latitude  is 
directly  proportional  to  the  vertical  gradient  of  the  zonal  angular 
momentum.  Since  9M0/3z  is  also  proportional  to  the  torque  exerted 
by  the  zonal  flow,  equation  (30)  may  be  interpreted  to  mean  that  the 
volume  transport  per  unit  time  in  the  meridional  direction  below  a 
certain  height  is  proportional  to  the  torque  exerted  by  the  zonal 
flow  above  that  height.  Since  there  is  no  volume  transport  below 
the  surface  z = 0,  the  zonal  flow  given  by  the  zero-order  equations 
cannot  exert  any  torque  on  the  earth’s  surface  in  the  steady  state. 

The  solutions  to  the  first-order  equations  can  be  obtained 
in  terms  of  a Green's  function,  but  as  these  solutions  cannot  be 
expressed  in  the  simple  form  found  from  the  zero-order  solutions, 
they  are  not  considered  in  detail  here.  However,  it  is  possible 
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to  drav  some  conclusions  about  the  torque  on  the  earth's  surface 
from  the  first-order  differential  equations.  Integrating  equation 
(16)  with,  respect  to  2 throughout  the  atmosphere  gives  the  result 
that  the  surface  torque  must  be 

T "■  - (1  - x2)  M 

*s  8x  x *0 

— z = 1 


(31) 


Since  M (l,x)  is  zero  at  the  noles  and  at  the  equator,  the  torque 
must  change  sign  in  the  interval  0 < |x  | <1.  This,  in  turn  implies 
that  the  zonal  flow  must  reverse  itself  in  some  region  between  the 
poles  and  the  equator.  As  this  reversal  of  the  zonal  flow  is  con- 
tributed by  the  first-order  equations,  and  consequently  is  small 
relative  to  the  zero-order  terms,  the  region  of  reversal  will  not 
extend  beyond  the  lower  layer  of  frictional  influence. 

8 • A Simple  Numerical  Example 

Most  of  the  features  of  this  model  can  be  illustrated  by 
considering  a numerical  example  in  which  dg0/8x  is  assumed  to  be  a 
simple  polynomial  in  x.  Since  gQ  is  functionally  dependent  on 
density  in  the  approximation  adopted  in  this  model,  assuming  8g0/8x 
to  be  known  is  equivalent  to  assuming  the  distribution  of  the  hori- 
zontal density  gradient  throughout  the  atmosphere  to  be  known.  The 
density  gradient  is  chosen  to  be  anti-symmetrical  about  the  equator, 
to  vanish  at  poles,  and  to  be  proportional  to 

= px3(l  - x2)  (32) 
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The  effects  of  friction  are  measured  by  tho  dimensionless  constant 
e,  which  in  this  example  is  chosen  to  h.vc  the  value  e = 10~3  to 
simplify  numerical  calculations. 

'’he  corresponding  zero-order  solutions  are 


*4>0  = - px(l  - x2)2  0(z,k) 

(33) 

M0  = px2(l  - x2)2  <X>(z,k) 

(34) 

gQ  = f(z)  + JL  x4(3  - 2x2) 

05 ) 

where  P is  a dimensionless  constant,  and  f(z)  is  an  arbitrary  func- 
tion of  z. 

Meridional  cross-sections  have  been  drawn  showing  isotachs 
of  the  zonal  velocity  component  in  figure  3,  and  streamlines  of 
the  meridional  circulation  in  figure  4.  The  normalization  factor 
for  the  zonal  velocity  component  is  umax  = 19.0  p m sec"*1,  and  for 
the  stream  function,  *P  = 6,85x10^ P ra^sec”1.  The  maximum  values  of 
the  velocity  components  obtained  in  the  numerical  example  arc* 

“max  ” 1?'°  0 ">  soc“1 

v __  = 0,59  P m sec"1  in  tho. upper  frictional  layor 

IQoa 

vmax  = P ® sec”1  in  the  lower  frictional  layer 

wmax  = P m sec-1. 

9 • Discussion 

The  simple  example  of  zonal  and  meridional  flow  given  above 
illustrates  the  type  of  flow  to  be  expected  from  this  model.  The 


flow  pattern  is  essentially  different  from  the  flew  to  to  expected 
or:  the  basis  ef  the  classical  Hadley  model  since  the  transport  of 
relative  angular  momentum  by  the  meridional  circulation  is  not 
taken  .into  account.  The  easterly  v/inds  near  the  equator  predicted 
by  the  model  are  not  due  to  horizontal  meridional  currents  as  sug- 
gested by  the  Hadley  model  but  arc  due  tc  a net  vertical  motion 
throughout  the  atmosphere  in  the  equatorial  regions.  The  easterly 
winds  have  an  appreciable  horizontal  extent  only  in  the  caso  that 

meridional  density  gradients  arc  very  small  in  the  equatorial  regions. 

2 

Thus,  a density  gradient  varying  as  x or  x in  the  vicinity  of  the 
equator  would  require  very  strong  east  winds  within  2 or  3 degrees 
of  the  equator.  Since  the  derivation  of  the  perturbation  equations 
implicitly  assumes  that  no  such  rapid  changes  in  velocity  occur, 
the  density  gradient  in  the  example  treated  above  is  chosen  to  vary 
as  near  the  equator. 

The  example  given  above  is  not  intended  to  resemble  the 
actual  circulation  in  the  atmosphere.  It  would  not  be  difficult 
to  choose  a density  distribution  which  would  give  a more  realistic 
zonal  velocity  profile;  however,  such  a choice  could  not  be  justi- 
fied from  the  dynamics  of  the  model  alone  and  would  require  the 
investigation  of  the  thermodynamic  processes  associated  with  the 
notion  for  its  justification, 

1 0 • Plans  for  Future  Work 

The  basic  purpose  of  pursuing  this  analysis  was  to  construct 
a model  in  which  heat  exchanges  could  be  considered  analytically. 


The  model  described  here  is  considered  as  the  first  approximation 
to  a more  realistic  model  and  would  have  to  be  general! cod  accord- 
ing to  the  problem  under  consideration.  Many  generalisations  can 
be  suggested;  for  example » it  may  be  necessary  to  take  into  account 
the  compressibility  of  the  atmosphere  before  significant  thermo- 
dynamic processes  can  be  considered.  Also  the  effects  of  lateral 
stresses  can  be  taken  into  account  by  considering  a nonisotropic 
eddy  viscosity.  Both  of  these  generalizations  make  the  mathematics 
of  the  model  more  complicated  but,  it  is  hoped,  not  impractical. 

The  effects  of  the  non-linear  terms  would  have  to  be  considered 
seriously  as  they  are  not  always  relatively  small  and  may  invalidate 
conclusions  made  from  the  linear  models.  It  is  not  proposed  to 
study  models  with  longitudinal  dependence  because  in  such  models 
time  dependent  terms  would  have  to  be  considered  and  the  analysis 
would  be  much  more  complicated. 

The  immediate  plan  is  to  apply  this  model  to  the  investiga- 
tion of  the  dynamics  and  thermodynamics  of  the  Antarctic  Circum- 
polar Current  because  only  minor  changes  of  the  present  model  are 
necessary  in  applying  it  to  the  ocean.  In  particular  it  would  be 
necessary  to  compute  the  solutions  to  the  zero-order  equations  cor- 
responding to  a given  distribution  of  surface  torque  exerted  by  the 
prevailing  winds. 
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Fig.  1.  The  function  t(z,k),  defined  bv  eq.  (2>) 
i s plotted  for  several  values  of  k.  For  the  numeric 
al  example  discussed,  this  function  is  proportional 
to  the  zonal  velocity  component  along  a vertical 
lino  throng!,  the  atmosphere. 


Flf.  2,  The  function  0(s,k),  defined  by  eq,  (26), 
!•  plotted  for  several  values  of  k.  For  the  numeric- 
al exanple  discussed,  this  function  la  proportional 
t*1*  stream  function  of  the  meridional  circulation 
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Fig.  3.  A meridional  cross-section  shoving  lsotacha  of  the  tonal  Telocit 
coup  orient  (eolld  lines)  normalised  to  unit  maxl.nura  speed  it  t • 1 The  sur* 
face  along  which  the  horltontal  pressure  gradient  vanishes  Is  marked  h»  « 
broken  line.  The  hatched  region  Indicates  the  hori»o-ital  extant  or  7 “ . 
erly  winds  as  given  by  the  first-order  equations.  extent  or  the  eaat 
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A meridional  cross-section  showing  streamlines  of  the  meridional 
i normalised  to  unit  maxisum  value  of  the  stream  function  at  * ■ 0.r. 
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